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We employ room-temperature ultrahigh vacuum scanning tunneling microscopy (UHV STM) and ab-initio 
calculations to study graphene flakes that were adsorbed onto the Si(lll)—7x7 surface. The characteristic 
7x7 reconstruction of this semiconductor substrate can be resolved through graphene at all scanning biases, 
thus indicating that the atomistic configuration of the semiconducting substrate is not altered upon graphene 
adsorption. Large-scale ab-initio calculations confirm these experimental observations and point to a lack of 
chemical bonding among interfacial graphene and silicon atoms. Our work provides insight into atomic-scale 
chemistry between graphene and highly-reactive surfaces, directing future passivation and chemical interaction 
work in graphene-based heterostructures. 


Graphene is a single layer of sp^-bonded carbon 
atoms^ that has an unusually large carrier mobil¬ 
ity and thermal conductivity.l^M^ Its electronic struc¬ 
ture is represented by a conical, relativistic-like energy- 
momentum relation at energies within 1 eV from its 
Fermi level. Silicon, on the other hand, remains the 
workhorse of present electronic technologies, forming 
a stable (111)—7 x 7 reconstructed surface after high 
temperature treatment Si (111)—2 x 1 surface 

reconstructiorP^®^ can only be produced by cleavage,^!^ 
giving way to the 7x7 reconstruction at 900 K. 

Many studies have addressed graphene’s electronic, 
mechanical, and topographic properties upon adsorp¬ 
tion onto semiconductin^i^HI^ and insulatin^^E^H^ sub¬ 
strates. Still, clean Si(lll) is perhaps the most reactive 
semiconductor to be interfaced with graphene or another 
inert two-dimensional material to date. It has a high 
density of dangling bonds, larger topographic modulation 
than Si(lOO)—2x1 or GaAs(llO), and higher crystalline 
ordering over Si02 or graphitized SiC. While the prop¬ 
erties of graphene adsorbed on Si(lll) were computa¬ 
tionally examined,!^ we aim to provide a comprehensive 
experimental and computational characterization of the 
intriguing graphene-Si(lll) system, thereby elucidating 
interfacial chemistry among graphene and reactive sur¬ 
faces. 

To this end, we interface the stable Si(lll)—7x7 
surface with graphene using the dry-contact transfer 
technique,!^ as illustrated in Fig. 1(a). Our transfer pro¬ 
cess produces a clean interface at the graphene-Si(lll) 
heterostructure, allowing examination with atomically 
resolved, ultrahigh vacuum (UHV) scanning tunneling 
microscopy (STM) and spectroscopy (STS). We frame 
these experimental results and explore the interfacial 
chemistry by ab-initio density-functional theory (DFT) 
calculations using a large supercell in the local density 
approximation (LDA). These calculations reproduce our 
experimental findings and reveal the absence of cova¬ 
lent G—Si bonds at the graphene-Si(lll) interface. Our 


article complements our previous STM studies of the 
atomic-scale interaction of carbon-based m ateria ls and 
technologically-relevant silicon substrates 

In the approach of Fig. 1(a), a fiberglass applicator 
collects graphene flakes and gently presses them onto a 
hydrogen-free Si(lll) substrate, allowing for subsequent, 
in situ, room-temperature STM interrogation. We note 
that the orientation of graphene flakes from the applica¬ 
tor is random, and, hence, we expect the registry between 
the graphene and Si(lll) lattices to be random also. A 
filled-states STM image of two graphene flakes on Si(lll) 
is displayed on Fig. 1(b); the graphene flakes are seen on a 
lighter color on the image. Despite the unknown relative 
orientation between graphene and Si(lll), the diamond¬ 
like structure characteristic of Si (111)—7 x 7 can be seen 
throughout the image. We superimpose two diamond- 
like features on the graphene and onto the bare Si(lll) 
surface to help make the area spanned by Si(lll) recon¬ 
structed supercells more evident. Such apparent registry 
regardless of relative orientation arises from the fact that 
the electronic density has a very large contribution from 
the dangling bonds on Si(lll)—even in scans acquired 
over graphene —in a n effect that we dubbed “electronic 
transparency. ”CSE1] 

We acquire the height profile of Fig. 1(c) along the line 
crossing the right side of the flake located on the left side 
of Fig. 1(b). It possesses an average height difference be¬ 
tween the graphene feature and Si(lll) of the order of 1.6 
A. Additionally, the height profile has an oscillatory pat¬ 
tern with a 27 A period, highlighted by vertical dashed 
lines. Si(lll) forms a triangular lattice with a lattice 
constant of 3.81 A. The 7x7 supercell has, in turn, a 
period of 7 x 3.81 = 26.7 A that matches the periodicity 
captured on the height profiles taken on the bare Si(lll) 
substrate and on top of the graphene flake. 

For a second graphene flake, we take an empty-states 
STM derivative image, shown in Fig. 1(d), and its corre¬ 
sponding height profile in Fig. 1(e). Both confirm the os¬ 
cillatory period of 27 A registered in Fig. 1(b). The aver- 
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FIG. 1. (a) Schematics of the dry contact transfer (DOT) technique, whereby a fiberglass applicator transfers graphene to a 
clean Si(lll)—7x7 surface in situ, (b) Filled states (—2 V) STM topograph of graphene flakes on the Si(lll)-7x7 surface; 
graphene flakes are seen on a whiter tone, (c) Height profile along the solid line on the STM topograph. The STM height of 
the graphene feature is ~ 1.6 ± 1.1 A. (d) Empty states (+1.9 V) STM derivative image for another graphene flake, (e) Height 
profile taken along solid line on subplot (d). Diamonds on the STM images indicate the area of a Si(lll)—7x7 supercell. 



FIG. 2. (a) Filled-states STM topograph and (b) derivative 
image of a graphene flake onSi(lll). The atomic arrangement 
of Si(lll)—7x7 can be clearly seen underneath graphene at 
the inset, (c) An empty states STM image and its derivative 
image help establishing that the 7x7 reconstruction of Si(lll) 
remains unchanged after graphene was adsorbed. 
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FIG. 3. (a) STM topograph in empty states of graphene on 
Si(lll). (b) The derivative image shows the relative misalign¬ 
ment among graphene and Si (111)—7x7 to be about 10^. (c) 
Another derivative image of a graphene flake on Si(lll)—7x7. 
(d) Average I — V data over the straight line in (c). No band 
gap is induced by the substrate or by the narrow width (17.1 
nm X 78 nm) of the flake. 


age STM height for this graphene feature is 1.7 A higher 
than the average height of the bare Si(lll) substrate, 
which is similar to the height difference estimated in 
Fig. 1(c), despite the different relative orientation among 
graphene and Si(lll) for these flakes. 

We now inspect subtle effects in the graphene- 
encapsulated Si(lll) substrate by bias-dependent scans. 
To this end. Fig. 2(a) shows a filled-states STM topo¬ 
graphic image of graphene on Si(lll), taken at a — 2 V 
bias and with a 100 pA tunneling current. Its derivative 
image is given in Fig. 2(b). These two images permit a 
clear identification of features characteristic of a Si(lll)- 


7x7 reconstructed surface. Furthermore, the 7x7 sur¬ 
face reconstruction is seen under empty-states scans, as 
confirmed in Fig. 2(c). The diamonds drawn on Fig. 2 in¬ 
dicate the unit cell of the silicon surface resolved through 
the graphene. The results from Fig. 2 help establish that 
the Si(lll)—7x7 surface is not altered chemically or mor¬ 
phologically upon graphene adsorption. 

In Fig. 3(a), we discern the relative orientation among 
graphene and the Si(lll)—7x7 surface by examining the 
graphene lattice over bright spots in STM images, where 
the STM resolves the graphene lattice rather than a con¬ 
volution with the substrate surface structure. There 
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is a relative misalignment of about 10*^ in the zigzag 
direction of the graphene lattice and the vertical line 
on the Si(lll)—7 x7 substrate for this particular flake 
(see derivative image, Fig. 3(b)). The yellow diamond 
superimposed in Fig. 3(a) reveals features characteris¬ 
tic of Si(lll), regardless of the incommensurate moire 
superlattice between the misaligned graphene and the 
Si(lll)-7 x7. Thus, even though graphene and Si(lll) 
may have incommensuration and thereby span a large su¬ 
percell, the Si(lll) STM image is not strongly convoluted 
by the moire superlattice. Consequently, theoretical cal¬ 
culations can reproduce features of the graphene-Si(lll) 
system with the smallest possible supercell, disregarding 
relative orientation concerns. This simplifying assump¬ 
tion from the experimental data in Figs. 1 to 3 makes 
comprehensive theoretical analysis possible, as detailed 
later. 

We show a derivative image of another graphene fea¬ 
ture in Fig. 3(c), and the straight line in the image corre¬ 
sponds to where we take averaged STS data. Fig. 3(d) 
gives these spectra, highlighting features consistent with 
graphene spectra obtained on other semiconducting sub¬ 
strates by STM.tl^®^ Any covalent interactions between 
the graphene and the Si (111)—7x7 substrate would dis¬ 
rupt the sp^ bonding network, leading to modifcations to 
the STS data which are absent here. 

It is time to discuss the computational study from 
Ref. |29j and a theory that explains our experimental 
observations. That work asserts that graphene can be 
placed on Si (111)—2 x 1 and compressed by more than 
11% to establish commensuration. In turn, unidirectional 
C—Si covalent bonds will be created in this extreme com¬ 
pression for graphene on Si (111)—2 x 1. These compu¬ 
tational results do not describe what we observe and are 
therefore inconsistent with our findings. 

Therefore, we follow our experimental observations 
which indicate that the features related to silicon dan¬ 
gling bonds are not hampered by a moire superlat¬ 
tice created at the graphene/Si(lll)—7 x 7 interface. 
Aiming for correctness and simplicity, we computation¬ 
ally examine the smallest graphene supercell matching 
Si(lll)-2 X 1. Fortunately, graphene has a lattice con¬ 
stant aog = 2.44 A, and a 11 x 11 graphene supercell 
matches the Si (111)—7 x 7 supercell within a negligible 
mismatch of 0.5 %. 

We create graphene/silicon supercells containing 687 
atoms, with a 10 A vacuum added along the vertical di¬ 
rection separating periodic images. De nsityT unctional 
theory calculations with th e SIES TA cod^^^I^ within the 
local density approximatiorP^^ yield reasonable separa¬ 
tions among weakly-bonded systems, giving reasonable 
results at the Si(l 11)/graphene interface. An equivalent 
mesh cutoff of 200 Ry is employed, and a conjugate- 
gradient ionic relaxation is performed until all force 
components became smaller than 0.02 eV/A. The re¬ 
sulting atomistic structure is shown in Fig. 4, and the 
graphene/Si(lll) supercell is emphasized by a yellow di¬ 
amond. The underlying Si(lll) reconstruction is quite 




FIG. 4. (a) Top and (b) side views of a supercell employed in 
calculations. The closest distance among silicon and carbon 
atoms is 2.1 A. The closer distance obtained experimentally 
(~ 1.7 A) may originate from the fact that the STM presses 
graphene down into the semiconducting substrate to establish 
the necessary feedback current. 



FIG. 5. (a) Filled and (b) empty-states iso-density surfaces for 
Si(lll)—7x7. Gomparison of these images to the experimen¬ 
tal STM images acquired underneath graphene in Figs. 2(b) 
and 2(c) helps confirming the lack of atomistic reconstruc¬ 
tion on the substrate upon the adsorption of graphene, (c) 
Empty-states iso-density surface for graphene on Si(lll). 


similar to the one originally obtained by Brommer and 
coworkers .1^ 

In our calculations, graphene adapts to the Si(lll) 
surface by bending vertically slightly. We note a verti¬ 
cal distance between graphene and the topmost Si atoms 
that varies between 2.1 and 3.1 A, for a 1.0 A variation 
in graphene height overall. The relative height among 
weakly bound materials can be further diminished un- 
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der the dynamical pressure exerted by the STM tipP^ 
thereby explaining smaller heights seen experimentally. 

Top views of the simulated density isosurfaces for 
Si (111)—7x7 are displayed in Figs. 5(a) and 5(b) and 
for the equilibrium structure containing graphene (c.f., 
Fig. 4) in Fig. 5(c). A clear superposition of graphene 
and Si(lll) features is apparent in Fig. 5(c). The elec¬ 
tronic hybridization among these systems increases when 
graphene is pushed in closer proximity to Si(lll), reveal¬ 
ing the semiconducting substrate’s features more clearly, 
consistent with our experimental results. Our computa¬ 
tional results show no covalent C—Si bonds between the 
graphene and Si(lll) in the equilibrium structure and, 
more importantly, they reproduce the periodicity seen in 
the experimental STM images. 

In summary, we studied monolayer graphene flakes ad¬ 
sorbed on the Si (111)—7x7 substrate with a combina¬ 
tion of STM measurements and ah-initio calculations. 
The characteristic 7x7 reconstruction of this semicon¬ 
ductor substrate was resolved through graphene in all of 
our scans, so its atomistic configuration was not altered 
upon the adsorption of graphene. No chemical bond¬ 
ing among graphene and silicon atoms occurred within 
this system. Therefore, we demonstrated that the typical 
7x7 reconstruction of Si(lll) remains unchanged after 
graphene adsorption on this highly reactive surface. Our 
report demonstrates graphene’s use as a protective layer 
for inherently reactive surfaces like Si(lll), facilitating 
ex-situ handling of these heterostructures. Additionally, 
these results help ongoing efforts in the integration of 
graphene with silicon for electronic and optoelectronic 
applications. 
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